Long-range proximity effects in superconductor-ferromagnet structures 
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We analyze the proximity effect in a superconduc- 
tor/ferromagnet (S/F) structure with a local inhomogeneity 
of the magnetization in the ferromagnet near the S/F inter- 
face. We demonstrate that not only the singlet but also the 
triplet component of the superconducting condensate is in- 
duced in the ferromagnet due to the proximity effect. Al- 
though the singlet component of the condensate penetrates 
into the ferromagnet over a short length ^ = yj D/h (h is 
the exchange field in the ferromagnet and D the diffusion co- 
efficient), the triplet component, being of the order of the 
singlet one at the S/F interface, penetrates over a long length 
U D/e (e is the energy). This long-range penetration leads to 
a significant increase of the ferromagnet conductance below 
the superconducting critical temperature T c . 

In recent experiments on S /F structures a considerable 
increase of the conductance below the superconducting 
critical temperature T c was observed Although in 

a recent work jlj it was suggested that such an increase 
may be due to scattering at the S/F interface, a careful 
measurement of the conductance demonstrated that the 
entire change of the conductance was due to an increase 
of the conductivity of the ferromagnet [Qj2). 

Such an increase would not be a great surprise if in- 
stead of the ferromagnet one had a normal metal N. It is 
well known (see for review J|,|||) that in S/N structures 
proximity effects can lead to a considerable increase of 
the conductance of the N wire provided its length does 
not exceed the phase breaking length L v . However in a 
S/F structure, if the superconducting pairing is singlet, 
the proximity effect is negligible at distances exceeding a 
much shorter length ~ This reduction of the proxim- 
ity effect due to the exchange field h of the ferromagnet 
is clear from the picture of Cooper pairs consisting of 
electrons with opposite spins. The proximity effect is 
not considerably affected by the exchange energy only 
if the latter is small h < T c . As concerns such strong 
ferromagnets as Fe or Co used in the experiments [EH, 
whose exchange energy h is by several orders of magni- 
tude larger than T c , a singlet pairing is impossible due 
to the strong difference in the energy dispersions for the 
two spin bands. At the same time, an arbitrary exchange 
field cannot destroy a triplet superconducting pairing be- 
cause the spins of the electrons forming Cooper pairs are 
already parallel. A possible role of the triplet compo- 
nent in transport properties of S/F structures has been 
noticed in Refs. @||, where the triplet component arose 



only as a result of mesoscopic fluctuations. However, in 
both cases the corrections to the conductance are much 
smaller than the observed ones. 

In this paper, we suggest a much more robust mecha- 
nism of formation of the triplet pairing in S/F structures, 
which is due to a local inhomogeneity of the magnetiza- 
tion M in the vicinity of the S/F interface. We show 
that the inhomogeneity generates a triplet component of 
the superconducting order parameter with an amplitude 
comparable with that of the singlet pairing. The penetra- 
tion length of the triplet component into the ferromagnet 
is equal to £ £ = y/ D/e, where the energy e is of the order 
of temperature T or the Thouless energy Et = D /L 2 , L 
is the sample size. The length £ £ is of the same order 
as that for the penetration of the superconducting pairs 
into a normal metal and therefore the increase of the con- 
ductance due to the proximity effect can be comparable 
with that in an S/N structure. 
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FIG. 1. Schematic view of the structure under considera- 
tion. In the inset is shown the structure, for which we calcu- 
late the conductance variation: two ferromagnetic wires con- 
nected to two ferromagnetic and two superconducting reser- 
voirs. 

We consider a structure shown in Fig.l and assume 
that the magnetization orientation varies linearly from 
a = at x — to a w — Qui at x = w. Here a is 
the angle between M and the z-axis. The case Qw = it 
corresponds to a domain wall with thickness w located 
at the S/F interface, while the model with the homo- 
geneous magnetization is recovered by putting Q = 0. 
This variation of M may also be brought about by an 
external magnetic field (see || and references therein). 
Of course, the variation of the magnetization considered 
is the simplest model of what may happen at the inter- 
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face in reality and we use it for simplicity. We consider 
the diffusive limit corresponding to a short mean free 
path. In this limit one may describe the S/F structure 
using the Usadel equation JTo| . The proximity effect in a 
S/F structure with an uniform magnetization M was an- 
alyzed in [jO]]. For the system with a non- homogeneous 
magnetization one should use a generalized form of this 
equation containing spin variables |L2"| |. The equation is 
non-linear and contains normal g and anomalous / qua- 
siclassical Green's functions. These functions are 4x4 
matrices in the Nambu<g>spin space. 

Assuming that the anomalous condensate function / 
is small, one can linearize the Usadel equation. This 
can be done if the transmission coefficient through the 
S/F interface is small due to a strong mismatch of the 
Fermi surfaces. Moreover, the order parameter in the 
superconductor can be strongly suppressed near the S /F 
interface if the transparency is high. In both cases one 
may assume a weak proximity effect, which presumably 
corresponds to the experiments. We write the Usadel 
equation for the matrix element (12) of / in the Nambu 
space. Then for the retarded matrix (in the spin space) 
Green's function f R we obtain (the index R is dropped) 



iDd 2 r f + 2ef - 2Aa 3 + (fV* + Vf 







(1) 



Here e is the energy, A is the superconducting order pa- 
rameter, which vanishes in the ferromagnet; &i are the 
Pauli matrices in the spin space, and the matrix V is de- 
fined as V = h (<r 3 cos a + (72 sin a) , where a varies with 
x as shown in Fig.l. This matrix describes the interac- 
tion between the exchange field and spins of the conduc- 
tion electrons and vanishes in the superconductor. Eq. 
(|l|) could be written for temperature anomalous Green 
functions f M at Matsubara frequencies u> , by replacing 
e — > i\u>\ and multiplying the last term by sgn(w) Jl2| . 
Eq. ([!]) is supplemented by the boundary conditions at 
the interface that can also be linearized |13|. Assuming 
that there are no spin-flip processes at the S/F interface, 
we have 



d x f = (p/R b ) fs , 

x=0 



(2) 



Here {...} is the anticommutator. In the region x > w, 
f n satisfies Eq. (||) with Q = 0. 

We see from Eq. (||) that the singlet component, com- 
muting with 173, and triplet component, anticommuting 
with 03, are mixed by the rotating exchange field h. In 
the region x > w the triplet and the singlet components 
decouple and their amplitudes should be found by match- 
ing the solutions at x = w. One should also use the 
boundary condition, Eq. (|^), and match the solutions 
in the ferromagnet and superconductor. It is clear that 
the singlet and triplet components of the anomalous func- 
tion /„ inevitably coexist in the ferromagnet. This fact is 
also known for the case of magnetic superconductors with 
Q / Q. In the region x > w the singlet part decays 
sharply but the triplet one survives over long distances. 
We are able to confirm these statements solving Eq. (||) 
with the boundary condition Eq. (|J). In the case of a 
homogeneous magnetization (Q — 0) the triplet pairing 
cannot be induced, which follows immediately from Eq. 
(||) connecting separately the singlet and triplet compo- 
nents at the opposite sides of the interface and Eq. (||). 

Eq. (||) can be solved exactly. The solution /„ can be 
written in the form 



/„ = a a A (x) + a 3 B (x) + ia x C (x) 



(4) 



The function C (x) in Eq. (|) is the amplitude of the 
triplet pairing, whereas the first and the second term de- 
scribe the singlet one. Substituting Eq. (|J) into Eq. (j^) 
we obtain a system of three equations for the functions 
A, B and C, which can be sought in the form 



A (x) = (Ai exp (-Kix) + Ai exp (k^x)) 



(5) 



The functions B(x) and C(x) can be written in a similar 
way. The eigenvalues Ki obey the algebraic equations 



-Q 2 )C -2(Qk)A = 



k*)B 



K 2 h A = 



(6) 



f) A - k\B + 2 (Qk) C = 



where p is the resistivity of the ferromagnet, i?b is the 
S/F interface resistance per unit area in the normal state, 
and f s = <7 3 A/Ve 2 - A 2 . The solution of Eq. (0) is 
trivial in the superconductor but needs some care in the 
ferromagnet. In the region < x < w the solution 
/ can be sought in the form / = U (x) f n U (x), where 
U (x) = a cos (Qx/2) +i<Ji sin (Qx/2). Substituting this 
expression into Eq. (0) and assuming that the solution 
depends on the coordinate x only we obtain the following 
equation for /„ 

- iDd 2 x J n +i (DQ 2 /2) (j n + aJn&^+DQ {dj n ,a^ 

+2ef n + h{a 3 J n }=0 (3) 



where n 2 = —2ie/D and n\ = —2ih/D (indices i were 
dropped). The eigenvalues k are the values at which the 
determinant of Eqs. (^) turns to zero. From the first 
equation of Eqs. (||) we see that in the homogeneous 
case (Q = 0) the triplet component has a characteristic 
penetration length ~ k^ 1 , but we see from Eq. that 
its amplitude is zero. If Q ^ 0, the triplet component 
C is coupled to the singlet component (A, B) induced in 
the ferromagnet according to the boundary condition Eq. 
(||) (proximity effect). If the width w is small, the triplet 
component changes only a little in the region (0, w) and 
spreads over a large distance of the order | 1 1 in the 
region (Q, L). In the case of a strong exchange field h, 
£h is very short <C w, £r), the singlet component 
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decays very fast over the length and its slowly varying 
part turns out to be small. In this case the first two 
eigenvalues Kip ~ (1 ± i)/£,h can be used everywhere 
in the ferromagnet (0 < x < L), where L is the length 
of the ferromagnet. As concerns the third eigenvalues, 
we obtain k 3 = -\/ K e + Q 2 in the interval (0,w), and 
«3 = k £ in the interval (w,L). The amplitude B% of the 
slowly varying part of the singlet component is equal to 
B 3 = 2 (QK 3 /K 2 h ) C 3 « C 3 . 

All the amplitudes should be chosen to satisfy the 
boundary conditions at x = (Eq. (|^)) and zero bound- 
ary condition at x — L. For the triplet component we 
obtain (we restore the indices R(A)) 

C R(A \x) = =R {QB{0) sinh (n e (L - x)) . 

[k £ cosh e cosh 83 + K3 sinh 9 e sinh 9 3 ] ^ (7) 

where w < x < L, B R ^ (0)=(p^ h /2R b ) ff (A) is the am- 
plitude of the singlet component at the S/F interface, 
Q c =k c L, Q^—k^w, and K^ A ^=y/^f2ie/D. One can see 
that the difference C R — C A is an even function of e. 
This is a direct consequence of the fact that C R — C A 
is proportional to the Fourier transform of the corre- 
lator K(t)=(ipi(t)ijj-i(0) + -0 T (O)^ T (i)>, which is even in 
time. In the Matsubara representation, C R ^ in Eq. (Q) 
should be replaced by C u with sgnw instead of (=f) and 
k^v/HAD, /s^^A/v 7 ^ 2 + A 2 instead of k? {A \ fg (A) 
respectively. Thus, C u corresponding to the temperature 
correlator K (r) = — < T T ip^ (0) ^-f (r) > is an odd func- 
tion of u> and the sum over all u> is zero in accordance 
with K (0)= K (0) =0. 

It is clear from Eq. (^), that the triplet component 
is of the same order of magnitude as the singlet one at 
the interface. Indeed, for the case w <C L we obtain 
from Eq. (fjj) |C(0)| ~ B(0)/ s'mha w , where a w = Qw 
is the angle characterizing the rotation of the magneti- 
zation. Therefore if the angle a w < 1 and the S/F inter- 
face transparency is not too small, the singlet and triplet 
components are not small. They are of the same order 
in the vicinity of the S/F interface, but while the singlet 
component decays abruptly over a short distance (~ 
the triplet one varies smoothly along the ferromagnet, 
turning to zero at the F reservoir. In Fig. 2 we plot the 
spatial dependence of the singlet |-B(x)| and the triplet 
|C(x)| components for two different Q. One can see that 
the singlet component decays abruptly undergoing the 
well known oscillations [fl5| while the triplet one decays 
to zero slowly. This decay in the region (0, w) increases 
with increasing Q. 

Thus, we come to a remarkable conclusion: the pen- 
etration of the superconducting condensate into a ferro- 
magnet may be similar to the penetration into a normal 
metal. The only difference is that, instead of the singlet 
component in the case of the normal metal, the triplet 
one penetrates into the ferromagnet. Of course, in order 



to induce the triplet component one needs an inhomo- 
geneity of the exchange field at the interface. 

The presence of the condensate function (triplet com- 
ponent) in the ferromagnet can lead to interesting long- 
range effects. One of them is a change of the conduc- 
tance of a ferromagnetic wire in a S/F structure (see 
inset in Fig.l) when the temperature is lowered below 
T c . This effect was observed first in S/N structures 
and later was successfully explained (see, e.g. reviews 
H[|). Now we consider the S/F structure shown in the 
inset of Fig.l. The normalized conductance variation 
SG = (G — G n ) /G n is given by the expression |l6[|: 



SG 



32T 



-Tr 



f R (x) - f A (x) 



(8) 



Here G n is the conductance in the normal state, Fy = 
1/2 [cosir 2 ((e + eV)/2T) + coslT 2 ((e - eV)/2T)] , and 
< .. > denotes the average over the length of the ferro- 
magnetic wire between the F reservoirs. The function / 
is given by the third term of Eq. (@) with C R = - (C A )* 
(we neglect the small singlet component). Substituting 
Eqs. (f| 0) into Eq. (|) one can determine the temper- 
ature dependence SG(T). Fig. 3 shows this dependence. 
We see that SG increases with decreasing temperature 
and saturates at T = 0. This monotonic behaviour of 5G 
contrasts with the so called reentrant behaviour of SG 
in S/N structures 0,^8| and is a result of broken time- 
reversal symmetry of the system under consideration. 

Available experimental data are still controversial. It 
has been established in a recent experiment || that the 
conductance of the ferromagnet does not change below 
T c and all changes in SG are due to changes of the S/F 
interface resistance R^. However, in other experiments 
Rb was negligibly small Q. The mechanism suggested 
in our work may explain the long-range effects observed 
in the experiments (!],§]. At the same time, the result 
of the experiment |J is not necessarily at odds with our 
findings. The inhomogeneity of the magnetic moment at 
the interface, which is the crucial ingredient of our theory, 
is not a phenomenon under control in these experiments. 
One can easily imagine that such inhomogeneity existed 
in the structures studied in Refs. |l],^] but was absent 
in those of Ref. ||]. The magnetic inhomogeneity near 
the interface may have different origins. Anyway, a more 
careful study of the possibility of a rotating magnetic 
moment should be performed to clarify this question. 

In order to explain the reentrant behaviour of SG(T) 
observed in Refs. [|]J|] one should take into account other 
mechanisms, as those analyzed in Refs. QJ^,|l^]. However, 
this question is beyond the scope of the present paper. 

We note that at the energies e of the order of Thouless 
energy e ~ Et the triplet component spreads over the full 
length L of the ferromagnetic wire (see Fig. 2). This long- 
range effect differs completely from the proximity effect 
in a ferromagnet with a uniform magnetization consid- 
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ered recently in Ref. P0[| . In the latter case the charac- 
teristic wave vector is equal to k± 2 = y/—2i(e ± h)/D (cf. 
Eqs. @). It was noted in Ref. @ that if e — > ±h, then 
^1,2 —> and the singlet component penetrates in the 
ferromagnet. If the characteristic energies e c h ~ Et,T 
are much less than h, the penetration length |/«i,2| 1 is 
of the order £/j and is much shorter then £t or L. 




FIG. 2. Spatial dependence of the singlet (dashed line) and 
triplet (solid line) components of | / 1 in the F wire for different 
values of a m . Here w = L/5, e — Et and h/Er = 400. 
Et = D/L 2 is the Thouless energy. 



0t w = 7T/10 



FIG. 3. The SG(T) dependence. Here 7 = p£, h /R b . 
A/E T > I and w/L = 0.05 

It is also interesting to note that a triplet component 
of the condensate function with the same symmetry (odd 
in frequency u> and even in momentum p) was suggested 
by Berenziskii ]2l]] as a possible phase in superfluid 3 He 
(this, so called "odd" superconductivity, was discussed 
in a subsequent paper |p2f ). Being symmetric in space, 
this component is not affected by potential impurities, 
in contrast to the case analyzed in Ref. f23|j , where the 
triplet component of the condensate was odd in space. 
While in 3 He this hypothetical condensate function is 
not realized (in 3 He it is odd in p but not in frequency), 
in our system this odd (in lu) triplet component does 
exist, although under special conditions described above. 

In conclusion, we have shown that in the presence of a 



local inhomogeneity of magnetization near the S /F inter- 
face, both the singlet and triplet components of the con- 
densate are created in the ferromagnet due to the prox- 
imity effect. The singlet component penetrates into the 
ferromagnet over a short length whereas the triplet 
component can spread over the full mesoscopic length 
of the ferromagnet. This long-range penetration of the 
triplet component should lead to a significant variation 
of the ferromagnet conductance below T c . 
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